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Abstract 

Large fluxes of low-energy cosmic rays in interstellar space 
may produce,, via unstable CNO beta-emitters, large fluxes of low- 
energy positrons and a detectable intensity of 0.51 MeV gamma rays. 
Even though on a galactic scale these cosmic-ray intensities are 
probably untenable since they conflict with the dynamics of the 
interstellar medium, they may exist on a smaller scale of the 
order of the stopping distances of both low-energy cosmic rays 
and positrons. We compare the results of a detailed calculation 
with direct positron measurements at the earth and we discuss the 
observability of the annihilation radiation above the isotropic 
x-ray background. 
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Low-Energy Cosmic Ray Positrons and 
0.51 MeV Gamma Rays from the Galaxy 


X. Introduction 

Measurements of the charge composition of cosmic-ray electrons 

12 3 

in the energy range of a few MeV to several GeV * * have indicated 

the existence of a detectable flux of cosmic-ray positrons which 

can best be interpreted as being of galactic origin. At lower ener- 

A’ 

gies, a measurement of Cline and Hones has placed an upper limit 

on the positron flux at about 1 MeV, with a strong possibility that 

a significant fraction of this flux is of cosmic origin. It has been 
5 

shown that the measured positron flux above about 100 MeV is con- 
sistent with that expected from the decay of positive pions produced 
in strong interactions between galactic cosmic-rays and the inter- 
stellar gas. 

Cosmic-ray positrons are produced not only from pion decay, 
but also from the beta-decay of unstable CNO nuclei produced pri- 
marily in relatively low-energy nuclear interactions. Since these 
interactions are nuclear transmutation, rather than pion-producing 
interactions, they involve relatively low threshold energies (of 
the order of 10 MeV rather than 300 MeV). Moreover, the beta-ray 
positrons are produced with a maximum energy of the order of 1 MeV, 
whereas the spectrum of pion-decay positrons peaks at about 35 MeV 
and decreases rapidly toward lower energies. Thus, if there exists 
a large flux of low-energy galactic cosmic-rays, these particles may 


3 



provide the most copious source of low-energy cosmic-ray positrons. 

This source then may contribute significantly to an observable flux 

near the earth at and below 1 MeV. Furthermore, because of their 

low energy, it is much more likely for beta-ray positrons to stop 

and annihilate in the galaxy to produce 0.51 MeV gamma-rays than 

the positrons from the pion decay process. 

Preliminary results on positron production from CNO beta- 

6 

emitters were reported by Ramaty and Misra and this problem was 

7 

considered in more detail by Misra . Estimates of positron pro- 

8 

duction from the CNO process were also made by Stecker , in con- 

9 

nection with galactic 0.51 MeV gamma-radiation, and by Verma , in 

an attempt to account for cosmic-ray electron measurements at low 
10 

energies . The formal treatment of the problem of galactic positron 

production in connection with galactic positron annihilation and the 

8 

resultant gamma-radiation was presented by Stecker . 

In the present paper we shall summarize the unpublished cal- 
7 

culations of Misra and relate them to the problem of 0.51 MeV gamma 
ray production. Furthermore, we shall discuss the energetics of the 
problem and relate th'*; results of our calculations to direct positron 
measurements at the earth and to the question of the observability of 
the annihilation radiation from the galaxy. 

II. Positron Production 

The production spectrum of radionuclei, q r (E^), per gram of 
interstellar material per second is given by 


4 



% N =«*Z ( £ i) n k ^K,Lk ( £ r; £ l) (1) 

i,k 

where and E r are the kinetic energies per nucleon of the pri- 
mary and secondary particles respectively, j is the intensity o^ 

i 

primary cosmic-ray nuclei of type i in interstellar space, n^ is 

the number density of interstellar gas atoms of type k, cr^ ^ is 

the interaction cross section for production of a radionucleus of 

type r and f (E *E ) is the normalized energy distribution function 
it 5 ik rl, i 

for production of a nucleus with energy E in an interaction of energy 


Among the variety of beta-emitting nuclei which may be formed 

XX 13 XA’ X5 

in galactic cosmic-ray interactions, only C , N ,0 and 0 con- 
tribute significantly to positron production. The principal re- 
actions and their threshold energies are given in Table 1 along with 

their decay modes and maximum positron energies. The cross sections 

11 

were summarized by Lingenfelter and Ramaty and by Audouze, Epherre 
12 

and Reeves . The distribution functions, f (E ;E.) are, in general, 

r , ik r * 1 

unknown. Their detailed form, however, is not important for the present 
problem since most of the positrons are emitted from nuclei decaying 
near rest. Thus, the positron production spectrum is determined mainly 
by the beta-ray spectrum of the parent nuclei. We have assumed that in 
these p-CNOinteractions, the kinetic energy per nucleon is equally 
distributed among all the secondary nuclei produced, so that 
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where A, and A, are the mass numbers of the incident and target 
i k 

nuclei respectively. 

We have assumed in these calculations that the ratios of 

H:C:N;0 in the interstellar gas are in the proportions 1:3,4 x 

10 ^:1.1 x 10’"^:5.8 x 10 based on the composition of the solar 
13 

atmosphere and which agree reasonably well with the universal 

14 15 

abundances gxven by Suess and Urey and Cameron 

The cosmic-ray energy spectrum and charge composition at solar 

minimum has been measured by a number of experimenters. The proton 

16 

spectrum above 20 MeV was summarized by Gloeckler and Jokipii and 

17 


below this energy is given by Fan et al. The CNO spectra above 

, 18 
50 MeV/nueleon were summarized by Meyer • Because of the lack of 

data at lower energies, we have assumed that these spectra are of 
0.4 

the form E^ , normalized to observations at higher energies. 

The interstellar primary cosmic-ray intensities can be obtained 
from the solar minimum spectra by multiplying by a modulation func- 
tion of the form 


_ rex? (!/«'?) 

/ ( e*/> *P) 

where R and |3 are rigidity and velocity respectively, is a 

characteristic transition rigidity that depends on the distribution 

of interplanetary magnetic field irregularities and 7] is a parameter 

which is space and time dependent and defines the total residual 

modulation. This modulation function was recently discussed by Ramaty 

19 

and Lingenfelter . They have shown that above a rigidity of 500 MV 


^ or R^ Ro 

R Ro 


(3) 
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2 3 

the observed cosmic-ray H' and He spectra are consistent with those 
obtained from nuclear interaction* in interstellar space, with “fl « 

2 

350 MV and a mean path length for galactic cosmic -rays of Xw 4 g/cm . 
Since there is no information on residual modulation at lower energies; 
we have to treat R Q as a free parameter. 

Since the lifetimes of the CNO beta emitters are short, the 
positron production spectrum, q^, can be obtained directly from the 
spectra q r , by taking into account the kinematics of the decay. Thus, 

where 7 r , 7 , and 7 are the Lorentz factors of the decaying isotope, 

t;Vi positron in the rest frame of the isotope, and the positron in the 

galactic rest frame, respectively; B is the maximum value of 7 and 

is given in Table 1 for the various isotopes involved. The function, 

P(7 ) is the normalized distribution function for positron production 

20 

from beta decay processes and is given by Fermi as 




The limits, 7“, are given by 
r 


- f J 1 (f -')*(?-$ 


n 



Using the data previously discussed, we have numerically 

solved equations (1) and (4) for the solar minimum spectrlN|and for 

T] = 350 MV with various values of R^. As an extreme assumption, we 

have also used a power-la w cosmic-ray spectrum of the form e' 2 ' 5 

with low-energy cutoffs, E^, at 5, 15 and 100 MeV/nucleon. 

The resultant production spectra together with the positrons 

5 

that would result from pion decay are shown in Figure 1. As can 
be seen, the CNO source becomes significant only if one assumes 
large fluxes of low-energy cosmic rays. The large positron pro- 
duction at low energies for = 200 MV and for the power law 

i ^ 14 

spectrum come* primarily from and 0 produced by low-energy 
protons via the reactions 0^ (p ,2p2n)N^ and N^(p,n)0^. 


III. The Galactic Positron Spectrum 

The positron density in interstellar space n^(^) can °^ ta ^ ne ^ 
from the production spectra given in Figure 1 by assuming a steady 
state and taking into account leakage from the galaxy, positron anni- 
hilation, and energy losses by ionization, bremsstrahlung, synchro- 
tron radiation and Compton collisions. We then solve the continuity 


equation 


A = j y,.) 




(7) 


where p is the density of interstellar material, r(^) t ^ ie tota l 

energy loss rate and j i s the effective positron survival time 

s 

against leakage and annihilation given by 


t 5 - 


( 8 ) 


# 
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The annihilation lifetime and various energy loss rates used 

8 

in this calculation were given by Stecker . For the low-energy 
positrons considered in this paper, ionization losses predominate 
so that «(f> is proportional to the gas density, p. 

The leakage lifetime can be expressed in terms of the mean 
amount of matter traversed by cosmic-rays, X 7 c ptf As mentioned 
above, X 4g/cm^. 

The solution of equation (7) for the equilibrium positron flux. 


J-j- - n^pc/A-p,- can be written as 


b»)- s ^ 


(9) 


where K(V) = r(j»)/ppe ^n cm^g 


is the stopping power of protons 


in hydrogen, and X is an equivalent annihilation-path- length given 

A 


by 



As can be seen, the positron intensity j + depends only on the 

interstellar cosmic ray intensity (through the production function 

and on the path length X. Furthermore, at low positron energy, where 

the range of these particles is much smaller than X, depends only 

on the cosmic-ray intensity itself. 

The positron fluxes j + were obtained by evaluating equation (9) $or 

2 

the production spectra given in Figure 1 and X = 4g/cm . The results, 

5 

combined with the calculations of Ramaty and Lingenfelter for positron 
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production via the rr-Ur-e process, are shown in Figure 2, together 
with the available positron measurements below 1 GeV. 

As can be seen, positrons with energies greater than about 
2 MeV come principally from pion decay. At lower energies the rela- 
tive contribution of the beta-emitters depends critically either on 
the value of (for ^ = 350 MV) or on the assumed low-energy cutoff 
(for a power law spectrum) . 

The energy densities of cosmic rays for the various spectra 

discussed above are given in Table 2, As can be seen from Figure 2 

and Table 2, unless the energy density in low energy cosmic rays is 

of the order of a few tens of ev/cm , most of the positron flux at 

~ 1 MeV would come from pion decay rather than CNO beta-decay and 

4 

would be small compared to the upper limit of Cline and Hones . If 

spread uniformly over the galactic disk, such energy densities lead 

21 

to serious difficulties regarding the stability of the galaxy . On 
the other hand, S since the ranges of both low-energy protons and posi- 

9 

trons are short ('-0.1 g/cm ), the ~ 1 MeV positron measurements sample 
only a small region of space which may not necessarily be representa- 
tive of the galaxy as a whole. 

For example, consider a 10 MeV proton; its range in hydrogen is 

2 3 
0.05g/cm which ,f or an ambient density of n atoms per cm leads to 

a lifetime t & 2 x 10"Vn years. The net streaming distance corre- 

1/2 

sponding to t is a « (4/3/1 |3Ct) , where rl is the mean free path 

1/2 

for diffusion. If ri = 1 pc, loo/n pc. Since the rate 

66 3 

of supernova explosions in the galaxy (volume « 2 x 10 cm ) is about 
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-2 1/2 
10 per year, in a spherical volume of radius 100/n pc a super - 

6 3/2 

nova would he expected to occur once every 2 x 10 n years. For 
3 

n :> 1 atom/ cm , this time interval is much larger than the lifetime 

against ionization given above. Therefore, if these protons are 

produced in supernova explosions and not by more frequent events such 

as novae or flare stars, their intensities would exhibit sharp maxima 

close to the time of the explosion and decay to much lower values 

22 

later on. According to Gold a supernova explosion may liberate 
52 

as much as 10 ergs. Assuming that 50% of this is in low-energy 

cosmic-rays, the mean energy density for about 2 x 10^7 n years over 

1/2 3/2 3 

a sphere of radius 100/n pc would be ~ 30 n er/cm . This is 

sufficient to produce a detectable flux of beta-ray positrons. On a 

galactic scale, however, the same sources of low energy cosmic rays 

-2 

occurring at a frequency of 10 per year would only produce an 
average energy density of 3/n ev/cm . For n = 1 this is somewhat 
large but not inconsistent with the overall energetics of the galaxy. 

Since the large fluxes of low-energy cosmic rays are restricted 
to small volumes and short times, they do not conflict with observa- 
tions that integrate over large distances, such as the temperature of 
HI clouds, based on 2f'cm observations, and the density of free electrons, 
based on the observed spectrum of non-thermal radio emission. (The 

connection of these quantities with low-energy cosmic rays was dis- 

23 

cussed by Balasubrahmanyan et al. ). The ^ 1 MeV positron measure- 
ments, however, by sampling only a small region of space corresponding 

o 5 1/2 

to their range of 0.2 g/cm , (t ps 10 / n years, a 70/n pc) may 
provide evidence for these locally enhanced fluxes of low-energy 
cosmic rays. 
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As can be seen from Figure 2, given a reasonable solar modula- 
tion, the positron measurements above a few MeV are all consistent 
with the calculated spectra from the rr-^-e process. Since the 

range of protons above pion production threshold is larger than 
2 

4g/cm , these spectra represent mean values over time periods com- 
parable to the positron leakage lifetime from the galaxy. The 
spatial and temporal inhomogeneities discussed above which would 
allow large but localized low-energy cosmic ray fluxes, do not apply 
to this energy domain. The positron spectrum, above a few MeV, ob- 
tained from a demodulated cosmic-ray distribution with T| = 350 MV, 
would therefore be a good representation of the interstellar positron 
intensity. 

IV. The 0.51 MeV Gamma -Ray Intensities 

The total number of positrons annihilating at or near rest per 
second per gram of interstellar gas is given by 



2 

Equation (11) was evaluated numerically for X = 4g/cm . The results 

are shown in Table 3 for the various spectra q^(^) given in Figure 1. 

As an extreme case, we have computed the annihilation rates that would 

result from an interstellar positron flux which is of the same order 

4 

as the upper limit given by Cline and Hones . Using a total positron 
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—2 2 

intensity of 2 x 10 particles of f MeV per cm.sec.sr we have com- 
puted a maximum annihilation rate corresponding to a positron life- 

5 

time against ionization of 10 /n years, and a mean rate corresponding 

6 3/2 

to an average time between supernova explosions of 5 x 10 (ft; years 

1/2 

in a spherical volume of radius 70/|t pc. These rates are also 
given in Table 3. 

g 

One of us (FWS) has previously shown that the positrons anni- 
hilate primarily from an S state of positronium with 75% of these 
annihilations producing a three-photon continuum rather than a two- 
photon line at 0.51 MeV. Therefore, on the average, one 0.51 MeV 
photon is produced for every two positrons which annihilate. Thus, 
the intensity of 0.51 MeV gamma rays observed along the line of 
sight as a function of galactic coordinates is given by 




( 12 ) 


where M(1 ,b " ) in g/cm is the amount of interstellar gas in the 

direction of observation. The resultant gamma-ray intensities for 

various directions of observation are given in Table 4. The values 
II II 

of M(1 ,b ') obtained from 21 cm observations are taken from Ginzburg 

24 II II 

and Syrovatskii . The values of M(1 ,b ) designated by "missing 

mass hypothesis" are based on the discussion of Stecker as ftwJ&i C 

26 

recent observations of 100 MeV gamma-rays from the galaxy as due to 
bremsstrahlung and tt° - production. 

8 II II 

As discussed in the previous paper > 1 0.5l (1 jb ) must bG 
greater than 3 x 10“^ cm “2 sr “l sec -l -j_ n order to be observable above 
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the X-ray background continuum by a detector with an energy resolu- 
tion of 5 keV which is the theoretical width of the 0 a f 51 MeV line. 

This intensity is therefore an absolute lower limit on the 0.51 MeV 
line intensity which must be present in order for the annihilation 
line to ever be detected. By comparing this lower limit with the 
calculated intensities given in Table 4, we see that the annihilation 
radiation could only be detected toward the galactic center, and, with 
the "massing mass hypothesis’’^ probably also as a disc average. The 
calculated gamma-ray fluxes, however, were obtained by assuming that 
the primary cosmic-ray intensities are spread uniformly along the 
line of sight over which the annihilation is formed. Such an assump- 
tion, for the demodulated cosmic ray spectrum with T| ~ 350 MV and 
R q = 200" MV and for the power law distribution, requires a mean cosmic- 
ray energy density of ~ 50 ev/cm^ in order to produce a detectable 
annihilation line. This energy density is at least an order of mag- 
nitude higher than that allowed by the general dynamics of the inter- 
21 

stellar medium . The energy requirements for the power-law spectrum 
with E„ ~ 100 MeV/nuc are smaller and therefore not necessarily 

Li 

inconsistent with the energy arguments mentioned above. A primary 

2 

cosmic-ray spectrum of this form, however, conflicts with the H and 

3 19 . 

He calculations and possibly with the positron measurements xn the 

3 

10 to 10 MeV region (see Figure 2), and therefore is probably not a 
good representation of the overall galactic cosmic-ray distribution. 

We conclude that the 0.51 MeV annihilation gamma-ray intensity, 
produced by positrons from both pion and CNO beta emitters for a 
homogeneous disc model will be smaller than the observed X-ray back- 
ground 27 and hence be unobservable. 
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As can be seen from Table 4, however, if the positron measure- 

4 

ment at 1 MeV of Cline and Hones is regarded as real flux rather 
than an upper limit and is spread uniformly along the line of 
sight, the annihilation line would be observable % both toward the 
galactic center and as a disc average. However, if these 1 MeV 
positrons are produced in supernova explosions and exhibit spatial 
and temporal inhomogeneities corresponding to their short ranges , 
the resultant gamma-ray flux would again be below the x-ray back- 
ground and unobservable. 

Since the galactic center is known to be an intense source of 

26 

high-energy gamma-rays and since the energy arguments that we have 
used do not necessarily hold for that region, the galactic center 
may be a detectable source of 0.51 MeV gamma rays. It may thus be 
more profitable to look for this radiation with a high-spatial- 
resolution detector than to look for a diffuse galactic flux. This 
argument is valid even if a diffuse flux is detectable, since as can 
be seen from Table 4, such a flux would still be more intense toward 
the galactic center. 
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Figure Captions 

1. Positron production spectra per gram of interstellar material using 
a power law cosmic ray spectrum with various low-energy cutoffs, E c , 
and demodulated solar minimum spectra with various values of and 

V 

2. Positron intensities in interstellar space together with the available 
positron measurements below 1 Bev. 
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Table 1 


Beta Emitter 
and Decay Mode 

Maximum Positron 
Energy (Mev) 

Production 

Mode 

Threshold Energy 
(Mev) 


0.97 

c 12 (p,pa) c 11 

20.2 

B 11 + 3 + + v 


N 14 (p,2p2n)C 1:L 

13.1 



0 16 (p , 3p3n)C 1:L 

28.6 

13 

N -> 

1.19 

N 1 ^ (p,p«l)N 13 

11.3 

c 13 + .p+ + v 


0^p,2p2n)N 13 

5.54 

t 

<}• 

H 

o 

1.86 

N^- 4 (p,n)0^ 4 

6 .4 

N 14 + p + + v 

i 



o 15 - 

1.73 

0 16 (p,pn)0 15 

16.54 

N 15 + P + + v 

i 

i 
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Cosmic Ray 


Solar 

Minimum 


l^= 350MV (E a 5 flucT) 
R 0 = 500MV R 0 = 200MV 


0.5 


0.6 


57 


Table 2 
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17 


50 







Table 3 


Positron Annihilation Rates (g**l sec"'*') 


Solar 

Minimum 

■m 

Power Law 

j + = 2xl0" 2 cm” 2 sec"*-sr~* - 

R 0 = 500MV 

R 0 = 200MV 

E 0 = X00 

Nucl 

n _ r Mev 
Ec “ 5 Nucl 

Mean 

Maximum 

l.lxio " 2 

1. 8xl0“ 3 

4. 3x10" 2 

3. 3x10” 2 

1.3X10” 1 

2. 2x10" 2 

5/2 

n 

0.5 
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